Massive Gauge Field Theory Without Higgs Mechanism 
III. Illustration of Unitarity 
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To illustrate the unitarity of the massive gauge field theory described in the foregoing papers, we 
calculate the scattering amplitudes up to the fourth order of perturbation by the optical theorem 
and the Landau-Cutkosky rule. In the calculations, it is shown that for a given process, if all the 
diagrams are taken into account, the contributions arising from the unphysical intermediate states 
included in the longitudinal part of the gauge boson propagator and in the ghost particle propagator 
are completely cancelled out with each other in the S-matrix elements. Therefore, the unitarity of 
the S-matrix is perfectly ensured. 
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, In the preceding paper called paper II, it was proved that the S-matrix given by the QCD with massive gluons is 
On ' independent of the gauge parameter. However, the gauge-independence of S-matrix previously was not considered 
' by some people to be a sufficient condition of the unitarity of a theory. Therefore, to demonstrate the unitarity, it 
is necessary to check whether the contributions arising from unphysical intermediate states to the S-matrix clement 
written for a given process are cancelled in a perturbative calculation . Historically, as mentioned in paper I^^l , 
several attempts'^! of establishing the massive gauge field theory without Higgs bosons were eventually negated. 
' The reason for this partly is due to that the theories were criticized to suffer from the difficulty of unitarity[^l~l^^l. 
T-H Whether our theory is unitary in perturbative calculations? That just is the question we want to answer in this paper. 

To exhibit the unitarity, we choose to calculate two-gauge boson and fermion-antifermion scattering amplitudes 
given in the perturbative approximation up to the order of g'^. According to the optical theoreml-'^^"^^! , we only need 
to compute the imaginary parts of these amplitudes. The imaginary part of an amplitude can be evaluated by the 
following formula'^^^^^l 



2ImTab ^Y^^acT^c (1-1) 

which was derived from the unitarity condition of S-matrix: S ~ S^S ~ 1 and the definition: S* — 1 + iT. We 
would like to emphasize that the above formula holds provided that the intermediate states {c} form a complete 
set. This means that when we use this formula to evaluate the imaginary part of an amplitude, we have to work in 
' Feynman gauge. In this gauge, the gauge boson propagator and the ghost particle one are given in the form 1^1 



and 



= .2 ^" 27- ' (1-2) 

tA-\k) - (1.3) 

where m denotes the gauge boson mass. 

In Eq. (1.2), the unit tensor g^j, can be represented by the completeness of the gauge boson intermediate states of 
polarization. 

3 

ff^. = I]e^(fc)eA.(fc) = P^..{k) + Q^,{k) (1.4) 

where P^v(k) and Q^u{k) are the transverse and longitudinal projectors, respectively. On the mass-shell, they are 
expressed as 

Ptj.u{k) = gf_,i, - kf_,k^/m'^, Qf^^{k) = k^^k^/n? (1.5) 
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It is noted here that in some previous works, I^^H^"''! the Landau gauge propagators were chosen at beginning to 
examine the unitarity through calculation of the imaginary part of transition amplitudes. This procedure, we think, 
is not reasonable and can not give a correct result in any case . This is because that in the Landau gauge, the 
gauge boson propagator only includes the transverse projector P^i, (k) which does not represent a complete set of the 
intermediate polarized states as seen from Eq. (1.4). Usually, the right hand side (RHS) of Eq. (1.1) is calculated by 
using the Landau- Cutkosky (L-C) rule'^^"^^!. By this rule, the intermediate propagators should be replaced by their 
imaginary parts 

Im{k^ - + ie)-i = -nd{kf - m^)0{ko) (L6) 

Utilizing the L-C rule to calculate the imaginary parts of the two-boson and fcrmion-antifcrmion scattering am- 
plitudes, we find, the unitarity of our theory is no problems. A key point to achieve this conclusion is how to deal 
with the loop diagram given by the gauge boson four-line vertex which was considered to give no contribution to the 
S-matrix clement in the previous investigations[''l[^^^^'''l. This diagram can be viewed as a limit of the loop diagram 
formed by the gauge boson three- lino vortices when one internal line in the latter loop is shrunk into a point. In this 
way, we are able to isolate from the former loop a term contributed from the unphysical intermediate states which 
just guarantees the cancellation of the unphysical amplitudes. 

The rest of this paper is arranged as follows. In section 2, we sketch the \mitarity of the S-matrix elements of order 
g"^. In section 3, we describe the calculations of the imaginary part of the two-gauge boson scattering amplitude in 
the perturbative approximation of the order g'^ and show how the unitarity is ensured. In section 4, the same thing 
will be done for the fermion-antifermion scattering. The last section serves to make comments and discussions. In 
Appendix, we will discuss the sign of imaginary parts of the loop diagrams by a rigorous calculation. 



2. UNITARITY OF THE TREE DIAGRAMS OF ORDER 

For tree diagrams of order g'^, the unitarity of their transition amplitudes is directly ensured by the on-mass shell 
condition. To illustrate this point, we discuss the fcrmion-antifcrmion (say, quark- antiquark) and two-gauge boson 
(say, two gluon ) scattering taking place in the S-channel as shown in Figs. (1) and (2). 

For the fermion-antifermion scattering, the S-matrix element may be written as 

Tfi = igY{Pi,P2)D^Ak)f{p'i,p'2) (2.1) 

where 

j''{pi,P2)=v{P2)Y^''u{pi) (2.2) 

and 

D,.{k) = " ^^1^.^' + "'"'H^' . ■ (2.3) 
fc^ — + te k-' — am'^ + le 

Noticing k = p\ + p2 = p'i+ p'2 and employing Dirac equation, it is easy to see 

V{PuP2) = (i (2.4) 

which holds due to that the theory only concerns the vector current and at each vertex the fermion and antifermion 

have the same mass as in the QCD. Therefore, the longitudinal term k^k^/k"^ in the propagator does not contribute 
to the S-matrix in the approximation of order g^ . In other words, the unphysical poles /c^ = and /c^ = arn? do not 
appear in the scattering amplitude. 

For the process depicted in Fig. (2), the transition amplitude is 



X r^,A (fci , k2 , q)D^^' {q)r^,,,x' ik[ ,k'^,q) 



(2.5) 



where 

'^ni'xikx, k2,q) = g^viki - k2)\ + gv\{k2 + g)^ - g\^l{kx + q)^ (2.6) 
and e'^{k) stands for the gauge boson wave function satisfying the transversity condition 
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(2.7) 



The transversity of the polarized states and the relation q = ki + k2 = k[ + k'2 directly lead to 

e''(fci)e"(fc2)r^.A(fci, fc2, q)q^ = 0. (2.8) 

This equality, analogous to Eq. (2.4), guarantees the removal of the unphysical poles from the S-matrix element 
written in Eq. (2.5). 

Similarly, for the t-channel and u-channel diagrams, it is easy to verify that the equalities in Eqs. (2.4) and (2.8) 
hold as well. These equalities ensure the S-matrix elements for these diagrams and other processes such as that a 
fermion and an antifcrmion annihilate into two bosons to be also unitary. 

The fact that the term k^ku/k"^ in the propagator gives no contribution to the S-matrix elements means that the 
S-matrix is gauge- independent at tree level. Therefore, in the gluon propagator, only the physical pole at k^ = rm? 
contributes to the S-matrix element and the gluon propagator given in the Fcynman gauge can reasonably be considered 
in calculation of the S-matrix elements. Certainly, the fact mentioned above allows us to write the intermediate 
states as transverse ones. When we evaluate the imaginary part of the transition amplitudes by the L-C rule, these 
intermediate states will be put on the mass shell. 



3. UNITARITY OF TWO-GLUON SCATTERING AMPLITUDE OF ORDER 



In the preceding section, it was shown that in the lowest order approximation of perturbation, the unitarity is 
no problem. How is it for higher order perturbative approximations? To answer this question, in this section, we 
investigate the unitarity of the two-gluon scattering amplitude given in the order of 5^. For this purpose, we only 
need to consider the diagrams shown in Figs. (3) and (4) and evaluate imaginary parts of the amplitudes of these 
diagrams. The diagrams involving fermion intermediate states are not necessarily taken into account because the 
fermion intermediate state is already physical. 

Fig. (3) contains eleven diagrams which have gauge boson intermediate states. Except for the last diagram shown 
in Fig. (3k), the other diagrams all have two-gauge boson intermediate states. If the unitarity condition is satisfied, 
in the amplitudes given by these diagrams, the unphysical parts arising from the longitudinally polarized intermediate 
states should be cancelled by the amplitudes of the five diagrams depicted in Fig. (4) which are of ghost intermediate 
states. To demonstrate this point, in the following, we separately calculate the imaginary parts of the amplitudes of 
all the diagrams in Figs. (3) and (4). 



A. The imaginary part of the diagrams in Figs.(3a)-(3j) 

By the L-C rule, the diagrams in Figs. (3a-3j) can be given by folding the tree diagrams shown in Fig. (5) with 
their conjugates. Through the folding, we obtain twice Figs. (3a-3f) and once Figs. (3g)-(3j). Noticing that each of 
the diagrams in Figs. (3g)-(3j) has a symmetry factor i, we can write the imaginary part of Figs. (3a)-(3j) as follows 

xT;;5^'^*(p;,p^;A:i,fc2)5'"''5""' ^ ' ' 



where 



and 



= ^^C^^Wp, +p2-k,- k2)nS{kl - rn?)e{k'i) , . 

X7r5(fc| - m^)0{kl) ^ ' ' 



T;t^''{pi,P2; fci, fc2) = E 4t'"'{pi,P2; fci, k2) 

r' (3-3) 
= -i5'e^(pi)e-(p2) E T^)^l"\pi,P2-M,k2) 



here T^v°''"^'^{pi,P2; ki, k2){i = 1, 2, 3, 4) stand for the matrix elements of Figs. (5a)-(5d) respectively. 
In light of Feynman rules, when we set 
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jace jbde jade jhce jabe jcde 



(3.4) 



where the superscripts abed of C, have been suppressed for simpUcity and notice the relations given by the energy- 
momentum conservation 



gi = Pi - All = fc2 - P2, q2 = Pi - k2 = ki - p2, qs = Pi + P2 = ki + k2, 
the functions Tpal^t'"^{pi,p2; fci, ^2) may be separately expressed as follows. 



For Fig. (5a), 



where 



and 



For Fig. (5b), 



where 



and 



For Fig. (5c), 



where 



and 



For Fig. (5d) 



where 



T^i^i^y{pi,P2;ki,k2) = ^p^+ferJ,^^A(Pi'^i'9i) 



xrii^^(P2,fc2,?i) 



Tp^aIpI'^I'^i) = ffpp(^i +Pi)x+ 9hx{Qi - ki)p - gxp{qi +Pi)^, 



^1-1^^2, k2, qi) = gav{P2 + k2)\ - Qvxiqi + k2)a + 5Aa(3l - P2)v 



T'^V^u"'^ {Pl , J)2 ; fcl , fc2 ) = ^T^^^T^^x (Pl ) ^2 , 92 ) X 



ri^j'^(P2,fci,gi) 



(2) 

^pAiPiM^ 92) = gpi^ipi + k2)\ + gux{q2 - k2)p - g\p{q2+Pi)u 



^alxiP'^' ^1' 92) = gp.a{P2 + ki)x - gpx{kl + q2)a + 5Aa(?2 - P2)n- 



T^aM J""' {pi,P2;ki,k2) = ^^plx ipi,P2,q3) 

xr^i'^(fci, A:2,g3) 



^faxiPl^P'i^ 93) = gpaiPl - P2)x + 5<ta(P2 + 93)p - gXp{q3 + Pl)a 

T^IaC^I' ^2, 93) = gnu{k2 - ki)x - g„x{k2 + qa)^ + gxpiqa + ki)^. 
Tl,tril''%i,P2; kuk2) = Ci7«^. + C2ifX^ + ^37^;^. 

TpCT/ni/ ~ gpcrgfiv ~ gtrngpvi 



^pJnu — gpcrgfiv gpfigavi 



(3.5) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 

(3.15) 

(3.16) 
(3.17) 



and 
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^p%u = 9p^9au - 9pu9ap ■ (3.18) 

The expressions in Eqs. (3.1)-(3.3), (3.6), (3.9) and (3.12), as indicated in Introduction, are all given in the Feyn- 
man gauge. When the intermediate states g^^ and g''" are decomposed into the physical and unphysical parts in 
accordance with Eqs. (1.4) and (1.5), Eq. (3.1) will be represented as 

2/mTi = i/ drT,?,5^'^T;^;^;^^*[P^^'(fci)P'^'''(fc2) + Q^^' {k^)9'"'' ,^ 

+,g^^'0'^'''(fc2)-Q^^'(fci)Q"'''(fc2)] ^' ' 

where P^^^' {ki) and Q^^' (ki) are defined in Eq. (1.5) and respectively represent the transverse (physical) and lon- 
gitudinal (unphysical) polarization intermediate states of gauge bosons. We see, except for the first term, the other 
terms in Eq. (3,19) arc all related to the unphysical intermediate states. These terms should be cancelled out in the 
total amplitude. In the following, we calculate these terms separately. In the calculations, we note, the transversity of 
the polarization vectors (see Eq. (2.7)), the relations written in Eq. (3.5) and the on shell property of the momenta 
Pi,P2, ki and k2 will be often used. 

1. Calculation of Tf,'"'Tp^^,'"^*Qi'i'' {ki)g""' 

To calculate the second term in Eq. (3.19), according to the definition in Eq. (3.3) for T^^'^''- and the expression in 

Eq. (1.5) for Q^'^'' {ki), we need to calculate the contractions k'^T^^Iit'"''{pi,P2;ki,k2){i = 1,2,3,4). Prom Eqs. (3.7) 
and (3.8). we find 

K^pixiPi' ki,qi) = -qipQix + 9Xp{ql - m') (3.20) 

and 

<li^^alx{P2M,qi) = -k2.qia- (3.21) 
Using these equalities, from Eq. (3.6), we obtain 

k'iTl'J^^''\pi,P2; fci, k2) = Cifc2.5W(gi) + C^V%{:p2, fe, ?i) (3.22) 

where 

gff(gi)= 2 ""'"l^- ■ (3-23) 
Similarly, from Eqs. (3.10) and (3.11), one can get 

^1 rSA(P2, fcl, 92) = -g2a52A + 9Xa{ql ' m^) (3.24) 

and 

92 rj,l(pi, fe, 92) = -q2pk2.. (3.25) 
Based on these equalities, it is found form Eq. (3.9) 

k'^T(l^-bcd^PuP2; kiM) = C2k2.SfX<l2) + C2r(l(pi, fc2, 92) (3.26) 

where 

SfJ{^2)= 2 . (3-27) 
Along the same line, we can derive from Eqs. (3.13) and (3.14) that 

K^^plxikiM, qs) = ki^qsx + kixk2. - 9.x{ql - m^) (3.28) 

and 
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(l3^fJx{PuP2, as) = (fci + k2)^r^p^M'P2, Qs) = 0, (3.29) 
thereby, we get from Eq. (3.12) 

k'^Tj,^-bc''{p,,P2; fei, fe) = C3fe.5(3)(pi,p2, fci) - C34iypi,P2, Qs) (3.30) 

where 

S^'J{pr,P2,h) = -^--l^fc^r2,(pi,P2,93). (3.31) 



In addition, from Eq. (3.15), we may write 

k'^T^a'i!^'!^'"^(pi,P2; ki, k2) = Ci^fX{ki) + C2lfJu{ki) ,^ 
+Ci^fl{ki) 

where 

ifaAki) = k^lfX^ = Qpcrkiu - Qp^kia, (3.33) 

ifJuiki) = A;i7p^U = 9pakiv - gcvkip (3.34) 

and 

ifXik^) = khfX'^ = 9<r.k^P - 9p.kia. (3.35) 
Summing up the results denoted in Eqs. (3.22), (3.26), (3.30) and (3.32) and noticing Eq. (3.3), we have 



fcKt(Pi,P2; fci, fc2) = -ig^e''{pi)e'^{p2)[k2.S'iy + G'ilf] (3.36) 



where 
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Sabcd ^ Y^CiS^^) (3.37) 



i=l 



and 



Gfcrf = Y.'^iG% (3.38) 



in which 



and 



G^'crl = ^%{P2, k2, qi) + li'Xiki), (3.39) 
= rfl{puk2, q2) + lfX{k2) (3.40) 

Gfl = -r(l(pi,P2, 93) + 7Sl(Pi,f>2, fci). (3.41) 



Employing the expressions given in Eqs. (3.8), (3.10), (3.13), (3.29) and (3.33)-(3.35) and considering the relations 
among the momenta as shown in Eq. (3.5) and the on-shell condition of the external momenta, it is not difScult to 
find the following relation 

G^pl - -G^pl = -Gfl. (3.42) 
Now, let us look at the color factors. According to the expression 

jabej^cde _ —(^Sachd " ^adhc) + {dacedbde " dbcedade) (3.43) 
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and defining 

Pi = ^SabScd + dabedcde, 

1^2 = JjSadSbc + dadedbce, (3.44) 
Pa — JjSacSbd + dacedbde, 

we may write 

Cl=Pl- P2, C2=Pl- P3, C3=/33- p2- (3.45) 

Substitution of Eq. (3.45) into Eq. (3.38) and use of Eq. (3.42) directly lead to 



+/33(4'l - 41) = 0. 



(3.46) 



This result makes Eq. (3.36) reduce to 
where 

gabcd ^ -,g^eP(p^)e''{p2)S;l^'^. (3.48) 
By using Eq. (3.47) and noticing fej = w?, we finally obtain 

T^t'^'^T^S'"^'' Q^"' {ki)g'"'' = S'^^^'^S"'^''"^* (3.49) 

It is emphasized that from the above derivation, we see, the four-line vertex diagram in Fig. (5d) plays an essential 
role to give the relation in Eq. (3.42) and hence to guarantee the cancellation of the second terms in Eq. (3.22), 
(3.26) and (3.30), which are free from the pole at = m^, as shown in Eq. (3.46). 

2. Calculation o/T^^''r;,'5"'''3''^'Q""'(fe) 

The procedure of calculating T"^t'^'^T'^,l,'^'^ g^^ Q""' {k2) completely parallels to that described in the former subsec- 
tion. From Eqs. (3.7) and (3.8), it follows that 

^2^^alxiP'^^ ^2, qi) = -qiaqix + QxMl - m^) (3.50) 



and 



qi'^i%iPuki,qi) = -h^qip- (3.51) 



These equalities allow us to get from Eq. (3.6) that 

^2^^ h,k2) = Ciki^Sl]) + CirW,(pi, fci, gi). (3.52) 

Based on the equalities 

f'2^?JxiPi^ ^2, q2) = -q2pq2X + Qpxiql - m?) (3.53) 

and 

92 r|5^(p2, fcl, 92) = -g2afcl/. (3.54) 

which are derived from Eqs. (3.10) and (3.11), it is found 

k^2T^%l""'{pi,P2; ki,k2) = C^ki^SfJ + C2T^^l{p2, k2, q2). (3.55) 
By making use of the equality 
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^2rjfjA(^i' ^2, qs) = -k2ixqi\ - kiixk2x + g^\{ql - m^) (3.56) 
which is derived from Eq. (3.14) and considering Eq. (3.29), we have 

fc2ri^ir.''''(pi,P2; Ai, fe) = Cafci^^^) + C3r(3)^(pi,P2, ga). (3.57) 
From Eq. (3.15), it is clear that 



+C2ji%{k2) + C3jill{k2) 

where 



v(2) 



(3.58) 



y^aU^2) = k^^l,l\, = gp„k2^ - g,^k2p, (3.59) 



Tp7,,{k2) = k^'y^^J^^ = gp^k2p - 5pA»^2a, (3.60) 
and 

7i%{k2) = fc27paU = 9pt^k2o- - gap.k2p- (3.61) 
Combining Eqs. (3.52), (3.55), (3.57) and (3.58), we obtain 

k^2Tf.fAPi,P2; kuk2) = -ig^eP(pi)e^ip2)[ki^Sf,''' + Gf,'^] (3.62) 

where Sfj"^ was defined in Eq. (3.37) and 

3 

Gfj^^{puP2; ki,k2) = Y.CiG% (3.63) 

i=l 

in which 

&pj, = r(;,Upi, fci, gi) + jl'J^{k2), (3.64) 

GfJ, = 4t (^^2, k2,q2) + l%{k2) (3.65) 

and 

Gf}^ = Tf}M^P2, qs) + l^X{k2). (3.66) 

Similar to Eq. (3.42), one may find 

&X = = -Gf},- (3.67) 
These relations and those given in Eq. (3.45) also lead Eq. (3.63) to vanish 

G^''>i.P2;fci,fc2) = 0. (3.68) 

Thus, Eq. (3.62) becomes 

fc2r;^''(pi,P2; fci, k2) = k.^s'^"-'' (3.69) 

where 5'"'"='^ was defined in Eq. (3.48) and thereby we have 

T^l'"^T^'^^,'"^*g'"^' Q""' {k2) = S^^^'^S"'^'"'^*. (3.70) 
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3. Calculation o/T^^''T;;^>'**Q'"''(fei)Q""'(A;2) 

To calculate T^t'"^TI^''t'''"^*Q^''' {ki)Q'"'' {k2), it is necessary to calculate k'^kl^Tfj"^. This may be done in several 
ways. For example, we may simply contract Eq. (3.69) with the vector k'( to give 

y(k^2Tfj'\puP2; fci, fe) = m^S"''^'^. (3.71) 

Certainly, paralleling to the procedure shown in the foregoing subsections, we may first compute k^ k^Tp^^f^"^ . For 
instance, by contracting Eq. (3.56) with fc^, wc derive 

kik^T%{ki,k2, 53) = -ki ■ k^qsx - m'^k^x + k^x{ql - m^). (3.72) 
Prom the above equality, noticing the identity in Eq. (3.29), it follows 

fcrfcaT^i^ir.''' = CWsfJ + C^k^,T%{pi.V2.qz). (3.73) 

The other terms can be given by contracting Eqs. (3.52), (3.55) and (3.58) with A:^. Summing all these terms, one 
can exactly obtain the result as written in Eq. (3.71). Employing Eq. (3.71), we get 

Tfj''^T^'y''''^*Q^'"' (ki)Q''''' (ki) = 3"'"='^ 8"'''''='^*. (3.74) 

Up to the present, the last three terms in Eq. (3.19) have been calculated. Inserting Eqs. (3.49), (3.70) and (3.74) 
into Eq. (3.19), we arrive at 

The second term in the above needs to be cancelled by the ghost diagrams. 

B. The imaginary part of the ghost diagrams in Fig. (4) 

The ghost diagrams in Fig. (4) can be given by folding the three tree diagrams plotted in Fig. (6) with their 
conjugates. The folding gives double Figs. (4a)-(4d) as well as Fig. (4e). Considering that the symmetry factor of 
Fig. (4e) is 1 other than i, the imaginary part of the transition amplitude of Fig. (4) may be represented as 

2lmTf-'^' = -y dTT-b'^{pi,p2;ki,k2)T-'''-d*{p[,p'^;ki,k2) 

- i / drT(3)«''-<^(pi ,p2;ku fc2)T(3)«'*''c* {p[ , p^; fci , ^2) ^ ' 

where 

3 

T'^'"'''{puP2;k,,k2) = 5^TW«''^''(pi,p2;fci,fc2) (3.77) 

i=l 

r(*)''*"^'^(pi,P2; All, fe) represent the matrix elements of Figs. (6a)-(6c), and the minus sign is inherent for the ghost 
loops. 

According to the Feynman rules and considering the transversity of the polarization states, it is clear that 

tW"^^''(Pi,P2; fc2) = S'W»'«='' (3.78) 

where 

Sii)abcd ^ -ig^eP{p,)e^p2)S^;} (3.79) 

here the sl^}{i = 1,2,3) were defined in Eqs. (3.23), (3.27) and (3.31) respectively. In accordance with Eq. (3.78), 
Eq. (3.76) can be expressed as 
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2ImTf"-'^' = J dTS^^'^^S"'^''"^* - rfT-5(3)a6cd^(3)a'6'cd* (3 gQ) 

When adding Eq. (3.80) to Eq.(3.75), we see, the second term in Eq. (3.75) is just cancelled by the first term in 
Eq. (3.80). However, still remains the second term in Eq. (3.80) which represents half of the contribution of the 
loop diagram in Fig. (4c) to the imaginary part of the amplitude, We are particularly interested in the fact that the 
first term in Eq. (3.80) contains the entire contributions from the ghost diagrams in Figs. (4a)-(4d) and half of the 
contribution of the diagram in Fig. (4e). They completely eliminate the unphysical part of the amplitudes given by 
Figs. (3a)-(3j), needless to introduce any extra scalar particle for this elimination. How to understand the remaining 
term in Eq. (3.80)? The occurrence of this term in the sum of the amplitudes given in Eqs. (3.75) and (3.80) is 
due to that the loop diagrams in Figs. (3g) and (4c) have different symmetry factors. Therefore, only half of Fig. 
(4e) is needed to cancel the unphysical part of Fig. (3g). It is reminded that until now, the loop diagram in Fig. 
(3k) has not been considered. This diagram, as Fig. (3g), has also a symmetry factor i and, as indicated soon later, 
gives a nonvanishing contribution to the scattering amplitude and its imaginary part in the case of massive gauge 
field theory. This contribution, of course, includes a part arising from the unphysical intermediate states when we 
calculate the loop diagram in Fig. (3k) in the Feynman gauge. Certainly, this unphysical part needs to be cancelled 
by the corresponding ghost diagram as well. From the theoretical logic, it is conceivable that the second term in Eq. 
(3.80) just serves such a cancellation. 



C. The imaginary part of the diagram in Fig. (3k) 

How to evaluate the imaginary part of the amplitude of Fig. (3k) by the L-C rule? This seems to be a difficult 
problem because we are not able to divide the diagram into two parts by cutting the internal boson line of the closed 
loop in Fig. (3k) without touching the vertex. However, we observe that when letting one boson line of the closed 
loop in Fig. (3g) shrink into a point. Fig. (3g) will be converted to Fig. (3k). This graphically intuitive observation 
suggests that the amplitude given by Fig. (3k) can be treated as a limit of the amplitude of Fig. (3g) when setting 
the momentum of one propagator in the loop shown in Fig. (3g) tend to infinity. In this way, we can isolate from 
Fig. (3k) the unphysical contribution which looks like to be given by two-particle intermediate states and hence is 
able to compare with the second term in Eq. (3.80). It is obvious that the difference between the both diagrams in 
Figs. (3k) and (3g) only lies in their loops, one of which is formed by the four-line vertex (See Fig. (7a)) and another 
by the three- line vertex (see Fig. (7b)). Therefore, it is only necessary to compare expressions of the two loops and 
establish a connection between them. 

The expression of the loop in Fig. (7a) is 

■n(X)°-bi N 2 racd rbcd\ i f 9^^ /0 01^ 

li-^y {q) = -g f f - gA^gA'.] J ^^^^4 _ ^ .^^ ■ (3.81) 

The imaginary part of the above function have been exactly calculated in Appendix. The result is 

imnir (g) = -^r'^f^^g,,, r '^^sinixm'). (3.82) 

V^T^j Jo ^ 

Clearly, it docs not vanish when the mass m is not equal to zero. The above result is derived in the Feynman gauge, 
therefore, contains the contribution arising from the unphysical intermediate states. 
The expression of Fig. (7b) will be written in the form 

^Til) = J + ^2 - 9)^'xl'\k„k2,q) (3.83) 

where 

^'-xl'\kl, k2, q) = igV''^'/'^'r^.A(fcl, fe, q) .3 

Xr^'.'v(fc2,fc2, (?)£""'' (fci)£>'^^'(fc2) ^ ■ ' 

in which the propagator D'^''{k) was given in Eq. (2.3) with the gauge parameter a = 1 and the vertex Tfj^^xiki, ^2, q) 
was defined in Eq. (2.6). Let us take the limit :|fc2^| — > 00. In this limit, the product of the propagator D"'^ (^2) and 
the vertices will approach to 
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-^9'"''[g^lvk2x - gxi,k2^j\[g^'u'k2\' - g\'u'k2,A (3.85) 

= --^[gtiti'k2\k2\' + 5AA'fc2/jfc2^' - g)j.\'k2\k2fj.' - gA^' ^'2m^'2A'] • 

If the tensor k2nk2v/kl behaves in such a way in the hmit 

k2nk2u/k2 Qnu (3.86) 

(This hmit will be justified in Appendix and a similar limitation for the polarization vector was proposed in the proof 
of 75-anomaly (see Ref. (20), Chapter 19)), then we find 

//-'^/^^^(.AV... - gV.v.) ,2_C:^,, (3.87) 



and hence 



^ p*") 

Particularly, in the physical region, the sign of the imaginary part of the amplitude n^^_^l"^((j) is the same as the 

corresponding part of the amplitude Il^^l^''{q), as will be demonstrated in Appendix. In view of the argument given 
above, the imaginary part of Fig. (3k) may equivalently be replaced by the imaginary part of Fig. (3g) in the limit 
\k2ij,\ 00. Thus, we can write 



t(i)«^ 



-i(3)abcdrp(S)a b ctZ* r 

+Q^'^'' {ki)g''"' + g^'^'' Q'^'^' {k2) - Q'^'''(A;i)Q"'^' (fca)]. 



2ImTt-'^' = ^JdT]im^,^^^^^4l)'''^X'l^^^^^^^ (3.89) 



The first term in the above only concerns the physical intermediate states. We do not pursue here what the limit for 
this term looks like because it is of no importance at present, we are interested in examining the other three terms. 
Look at the expression given in Eq. (3.28). The first term in it can be ignored due to the equality in Eq. (3.29). 
The last term can also be neglected comparing to the second term in the limit |fc2/i| — > oo. Thus, Eq. (3.30) will be 
reduced to 

fcfl^j?;.''' - Csk2.SlS (3.90) 

where Spa was defined in Eq. (3.31) and is irrelevant to k2. The result in Eq. (3.90) enables us to write the second 
term in Eq. (3.89) as 

where S^^^'^'"^'^ was defined in Eq. (3.79). In the above, the compatibility of the on-shell condition fc^ = with the 
limit |fc2p| 00 has been noticed . 

By the same reason as stated above, only the second term in Eq. (3.56) should be considered in the limit. Therefore, 
Eq. (3.57) is approximated to 



where 



|fe2M|->oo ^3 — + IS 

With the result in Eq. (3.92), the third term in Eq. (3.89) becomes 



k^2T^%l""' - -C^k^y^Sf) (3.92) 



5(3) = lim ^"^P^]:^^^^. (3.93) 



|'S2^|^00 

where 
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Si3)abcd ^ _ig2e^(p,)e-(p2)C35(3). (3.95) 

Similarly, in the limit \k2fi\ oo, we can neglect the first term (due to Eq. (3.29)) and the last term in Eq. (3.72). 
The second term in Eq. (3.72) permits us to rewrite Eq. (3.73) in the form 

Kk^^T^lll'^" - -CWsfJ (3.96) 
which may more directly be given by contracting Eq. (3.92) with fcf. From this result, it is clear to see 

lim Tl^l^''^'^Tl^Jlf^''''^*Qf'f'\ki)Q''''\k2) =S^^^^^^^ (3.97) 

|fe2M|-»°0 

On inserting Eqs. (3.91), (3.94) and (3.97) into Eq. (3.89), we see, the last two terms in Eq. (3.89) are cancelled 
with each other. As a result, we have 

2ImTi''-''' = i / dTlim|fe,,^|^^ri3.)''^^'^T^?)?''''^''*P'"'' {k,)P-^-^' {k^) 

j^l j ^^gCi)ahcdg(?,)a'b'cd* l,0.yo; 

The second term above is just cancelled by the second term in Eq. (3.80). Combining the results given in Eqs. 
(3. 75), (3. 80) and (3.98), we obtain the total amplitude as follows 

3 

2j^ya&a b ^ j2 27mTf^« = 1 / dTT^l^'^T^^.y'^* P^'^' {ki)P'"' (^2) 

+i / rfTlim|fe,^|^^T^e^"''^'^T^?)?''''^''*P'"'' {kr)P'^-' (fc2) 

which is only related to the physical intermediate states. Thus, the proof of the unitarity is accomplished. 

We note here that the results given in this subsection rely on how to correctly treat the limit procedure. As will be 
shown in Appendix, the limit given in Eq. (3.86) is the only choice of converting Fig. (3g) into Fig. (3k) when the 
relation in Eq. (3.29) is considered. Similarly, to obtain the desirable limiting results presented in Eqs. (3.91), (3.94) 
and (3.97) , the reasonable expressions in Eqs. (3.28), (3.56) and (3.72) are necessary to be used. 



4. UNITARITY OF FERMION-ANTIFERMION SCATTERING AMPLITUDE OF ORDER 

In this section, to illustrate the unitarity of the theory further, we evaluate the imaginary part of the fermion- 
antifermion scattering amplitude in the perturbative approximation of order g'^. For this purpose, it is only necessary 
to consider the diagrams shown in Fig. (8). 

The diagrams in Figs. (8a)-(8e) can be reconstructed by folding the tree diagrams in Figs. (9a)-(9c) with their 
conjugates. Since the folding gives two times of Figs. (8a)-(8d) and one time of Fig. (8e) which possesses a symmetry 
factor i, the imaginary parts of the amplitudes given by Figs. (8a)- (8e) may be represented as 

= I / drT^lT^}':, [P''"' {k^)P'^-' (fc2) + Q""' (fcOff-^-' (4.1) 

where 

3 

T;t{Pi,P2;ki,k2) = YTl:}''\pi,P2-M,k2) (4.2) 

i=l 

T^v'^''{pi,P2; fei, fe) {i = 1, 2, 3) denote the matrix elements of Figs. (9a)-(9c) respectively. According to the Feynman 
rules, they can be written as 

\b b- _l_ Ji/T \*i 

T(,l^''\p„P2;k„k2) = -^g^^p2)-J^ ^^^ - M^ + ie Y^^''^P'^^ ^^"^^ 

\ a u n -I- A/f \^ 

Ti^r\pi,P2; h,k2) = -ig\{p2)-^, -P2r-M2 + is T>^(p^) (4.4) 
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and 

T^l^''\pi,P2;h,k2) = - ,2^- r^.x(ki,k2,q)vip2)^lMpi) (4.5) 
^ — m'^ + ie 2 

where T ^^\{ki, k2, q) was defined in Eq. (2.6), M is the fermion mass and p =p^^/ii. 

For evaluating the second term in Eq. (4.1), we need to compute the contraction of T^''}"''' with k^. By applying 
Dirac equation, the on-mass shell condition and the relation q = ki + k2 = Pi + P2, one may get 



= g^r^^v{p2)^^MPi) 



2 ' 2 J "^^yy^j (4.6) 
= y~ j'""v\:p2)-y1v''j-\pi) 

and 



where 



KTlS^' = -9'r'"v{P2)^lMPi) + fc2.5«^ (4.7) 



~2 £abc \c 

= \ 2 ^(P2)^ki»(pi) (4.. 

2pi • P2 + 2 



Adding Eq. (4.7) to Eq. (4.6), we find 



fc^T;,^ = fc2.5«^ (4.9) 



As seen from the above, there is a cancellation among the diagrams in Figs. (9a)-(9c). From Eq. (4.9), one may 
derive 

T^IXtl'^' Q^^' (fci ) g""' = S'^^S''^* . (4.10) 
Let us turn to calculate the third term in Eq. (4.1). Along the same line stated above, one may get 

k^lTj^)'^' + T(2)«f] = -g^p»c^p,)^^^u{p^) (4.11) 



and 



where 



^2^(3;"" = 5'r'MP2)y7M«(Pi) + h^S'^' (4.12) 



= 2 " '2 , . V{P2)^^2U{P^). (4.13) 



q-' — m-' +ie 2 



Prom the equality 



it follows that 



^^(?52)y(ki + k2)w(pi) = t'(p2)y(Pi +P2)u(pi) = 0, (4.14) 



S"'' = S"''. (4.15) 
Adding Eq. (4.12) to Eq. (4.11) and noticing Eq. (4.15), we have 

k2T;t = fci^S"'. (4.16) 

This result gives rise to 

T;^T;,f ' 5""' Q"""' (fc2 ) = S^^S"^* (4.17) 
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For evaluating the last term in Eq. (4.1), we may use the following equalities which are obtained by contracting 
Eqs. (4.11) and (4.12) with k'^ 

Kk^lT^l^-" + T(2)«f] = -g^r'^vip^) ^i^MPi) (4.18) 



and 



Kk2TlS)^' = 5'r'MP2)ykiw(pi) + m'^S'^K (4.19) 



These equalities and Eq. (4.15) lead to 



fcf fc^^rt = m^S"''. (4.20) 



This result allows us to write the last term in Eq. (4.1) in the form 

TlitTlitQ^^' {ki)Q'"'' {k2) = S'^^S''^*. (4.21) 
Substituting Eqs. (4.10), (4.17) and (4.21) in Eq. (4.1), we arrive at 

2/mTi = ^J drT^^^T^f^.P^'^' {ki)P'"'' {k2) + dTS"''S"^*. (4.22) 

The ghost diagram in Fig. (8f) can be given by folding the tree diagram in Fig. (9d) with its conjugate. Therefore, 
the imaginary part of Fig. (8f), by the Feynman rules, can be written as 

2ImT2 = - j dTS'^^S"''* (4.23) 

In complete analogy with the two-boson scattering discussed in the preceding section, the second term in Eq. (4.22) 
can only cancel half of the above amplitude. The reason for this still is due to the difference between the symmetry 
factors of Figs. (8e) and (8f). To achieve a complete cancellation, it is necessary to consider the contribution of the 
diagram in Fig. (8g) . This diagram can also be treated as a limit of the diagram in Fig. (8e) when the momentum 
of one internal line in the loop tends to infinity, 

2ImT, = i/drlim|fe,^|^^Tg)"''r(?)f''*[P'^'^'(fci)^''''''(fc2) .4 24) 

+Q"^'(fci),9""' +,9^^'Q""'(fc2) - Q^^'(fci)0""'(fc2)]. 

In the limit:|fc2/i| — > 00, comparing to the second terms in Eqs. (4.7), (4.12) and (4.19), the first terms in these 
equations can be ignored. Thus, Eqs. (4.7), (4.12) and (4.19) respectively reduce to 

^/.j<3)af. ^ k2uS''^. (4.25) 

k^Tj^J'^' « (4.26) 

and 

k'^k^Tl^J"^ « M'^S"^ (4.27) 
On inserting these expressions into Eq. (4.24), we have 

2ImT, = \JdT)iu,^,,^\^^TlS)'''Tpf*P^^\k,)P^^\k2) .428) 

Thus, as shown before, it is indeed possible to find a way which allows us to isolate from Fig. (8g) the unphysical 

part of the amplitude like the second term in Eq. (4.28). 

Summing the results denoted in Eqs .(4.22), (4.23) and (4.28), we finaly obtain the imaginary part of the total 
amplitude such that 

2ImT = \ J dTT-lT;f*P^^^^' {kr)P--' {k2) 
+\!dTY.ra\,^^\^^TlS'''T^]:^'*P'^^'{k,)P^^'{k2) ^ ' 

in which the unphysical contributions are all cancelled. Thus, the unitarity is ensured. 
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5. COMMENTS AND DISCUSSIONS 



In the previous sections, the unitarity of our theory has been illustrated by evaluating the imaginary parts of 
two-gauge boson and fcrmion-antifcrmion scattering amplitudes up to the fourth order perturbation. The imaginary 
parts of the amplitudes were calculated by means of the L-C rule. In this kind of calculation, we have to first 
work in the Feynman gauge because the formula used requires the intermediate states to be complete. The gauge 
boson propagator given in the Feynman gauge contains unphysical longitudinal intermediate excitations though, it 
has been proved that these unphysical intermediate states are eventually cancelled in the S-matrix elements, leaving 
only physical transverse intermediate states in the S-matrix elements as given by the gauge boson propagator written 
in the unitary gauge 

D^.{k) = ^rr^^^^""' - (5-1) 

Although the unitarity of the S-matrix elements is proved in the Feynman gauge, it would be true for other gauges 
since it has been exactly proved that the S-matrix is gauge-independent. 

As mentioned in Introduction, In the previous works of examining the unitarity of some kinds of massive non-Abelian 
gauge field theoriesI^'^H^^l, the gauge boson propagator in the Landau gauge 

was chosen to calculate the imaginary part of scattering amplitudes by the L-C rule. For such a calculation, as 
mentioned in Introduction, the Landau gauge actually is not suitable since the intermediate states characterized by 
the transverse projector appearing in the numerator of the above propagator does not form a complete set. The 
unsuitability of the procedure may be seen from the massless gauge theory. The unitarity of the theory was exampled 
by c;omputing the imaginary part of the fermion-antifermion scattering amplitude of order (/^ in the Feynman gauge'^^l. 
However, if one tries to perform the proof in the Landau gauge, he could not get a reasonable result. Particularly, 
the longitudinal projector fc^fc^/fc^ in Eq. (5.2) where m = could not be given an unambiguous definition on the 
mass shell fc^ = since the momentum k on the mass shell becomes an isotropic vector. 

Another point we would like to emphasize is that for examining the unitarity of a massive gauge field theory, it is 
only necessary to evaluate the S-matrix element between the physical transversely polarized states of gauge bosons. 
In this way, it was shown in sections (2) and (3) that the unitarity is well satisfied. Particularly, the calculation in 
section 3 indicates that in all the diagrams depicted in Figs. (3) and (4) except for the loop diagrams involving Fig. 
(3k) and a part of Fig. (4e), there is a natural cancellation among the contributions coming from the unphysical gauge 
boson and ghost particle intermediate states, without the help of any scalar particle. This result and the theoretical 
logic strongly suggest that the same cancellation between the unphysical contributions arising from Fig. (3k) and a 
part of Fig. (4c) is definite to happen. To achieve this cancellation, the loop diagram formed by the gauge boson 
four-line vertex is necessary to be recast in the form as if it is given by the two-particle intermediate states so as to 
be able to compare its contribution with that given by the other loop diagrams. For this purpose, we proposed in 
section 3 a reasonable limiting procedure which allows us to reach the cancellation mentioned above. The results we 
obtained are undoubtedly correct. It should be noted that in all the previous investigationsI^^^I^^^I^^l of the unitarity 
problem, the loop diagrams given by the gauge boson four-line vertex such as Figs. (3k) and (8g) were never taken 
into account in the cancellation of the unphysical amplitudes. From the calculations described in sections (3) and 
(4), it is clearly seen that these loop diagrams play an essential role to guarantee the cancellation of the unphysical 
part of the amplitudes and hence the unitarity of the S-matrix elements. It is mentioned that the theories presented 
in Refs. (2) and (6) which were pointed out to be non-unitary by Mohapatra et. al. and some othcrsl^^^^-^^l are not 
correct because the Feynman rule concerning the closed ghost loop has an extra factor i other than 1 as given in 
our theory. In this paper, the unitarity has been proved in the perturbation approximation up to the order of g^. 
For higher order approximations, we believe that for a given process, if all the diagrams are taken into account and 
treated appropriately, the unitarity would be proved to be no problem. 

It should be noted that in this paper and the former papers, we only limit ourself to discuss the theory ( for example, 
the QCD with massive gluons) in which all the gauge bosons are required to have the same masses. This requirement 
is necessary to make the theory to be gauge- invariant and unitary just as the case we met in nuclear physics where the 
nucleon-pion interacting system is of SU(2)-symmetry provided that the masses of all the pions are considered to be 
the same and the mass difference between proton and neutron is ignored. .Similarly, for the weak interaction theory in 
which apart from the vector currents of fermions, the axial vector currents of fermions are included as well, certainly, 
we may build up a SU(2)- symmetric theory without introducing the Higgs mechanism. The action of this theory is 
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constrained by the Lorcntz condition and may contain the gauge boson mass term in it. But, the gauge-invariance 
of the theory requires that all the gauge bosons must be of the same mass and the fcrmions must be massless. In 
this case, the unitarity of the theory, as easily checked, is no problem. Nevertheless, if the charged and neutral gauge 
bosons are required to have different masses and charged fermions are massive, it is inevitable to work with the theory 
constructed by means of the Higgs mechanism. 
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7. APPENDIX: EXAMINATION OF SIGNS OF THE IMAGINARY PARTS OF THE LOOP DIAGRAMS 

It was mentioned in section (3) that the imaginary part of the matrix clement of the loop in Fig. (7a) has the same 
sign as that given by the loop in Fig. (7b) at least for the large momentum k2. To convince ourself of this point, we 
investigate the imaginary part of the loops in a parametric representation. In this representation, the propagator will 
be expressed in the form ^^''1 



1 

-i 



daie'"'('=^-'"'+'^). (Al) 







kf — + is 

With this representation, the matrix element shown in Eq. (3.81) for Fig. (7a) may be rewritten as 

^xx'^iQ) = -i9^r'"'f'""'gxx'J^'^ (A2) 

where 

J(l) = J .d^k 1 



^r^^/TiiT'^'"^' (A3) 



<fk — ■ 
_ 1 r°° da ^-ia(m -is) 

In the above, we have used the representation in Eq. (Al) and the formula of Fresnel integral!^''! 



I 



e'"*^ = — (A4) 

{2n)^ (47ra)2 ^ ' 



For the loop in Fig. (7b), we confine ourselves to investigate its expression in the limit of large momentum k2 for 
the purpose of comparing to the one given in Eqs. (A2) and (A3). As stated in Eqs. (3.85)-(3.88), in order to convert 
the matrix element of Fig. (7b) to the one for Fig. (7a), it is necessary to take the approximate expression denoted in 
Eq. (3.85) and the Hmit assumed in Eq. (3.86) which is applied to Eq. (3.85). To achieve such a conversion, as easily 
verified, instead of Eq. (3.85), we may simply take an equivalent approximate expression such that 

g^^'g'"''T^^x{ki,k2,q)T^,^,x'{ki,k2,q) ^ 6k2xk2X'. (A5) 
With this expression, Eqs. (3.83) and (3.84) may be rewritten as 

n[if{q) « 2,g^r'f''jfl (q) (A6) 

where 

J^xlil) = I Wf^wf^i^-m^ - /Cl - k2) 

Employing the parametrization given in Eq. (Al) and the Fourier representation of the (5-function, Eq. (A7) reads 

•^fizfiz f d^fe „»«2fc2-|f--i(a;i+a2)(TO^-fe) 
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Upon completing the integrations over fci, k2 and z by using the formula of Presnel integral, we have 

t(^) _ i r°° daida2 \„^^^ i _ 1 1 



where 



Inserting the identity 



/■oo 

/ c?a;(5(a; - ai - = 1 (All) 
Jo 

into Eq. (A9) and then making the transformation a, — > araj, one can write 

Noticing the equalities shown in Eqs. (3.29) and (4.14), the second term in the parenthesis, actually, can be ignored in 

(2) 

the scattering amplitude. Thus, the function J)^^ (q) is only proportional to the unit tensor gw . This result precisely 
justifies the limit taken in Eq. (3.86). On substituting Eq. (A12) into Eq. (A6), and performing the integration over 
a2, one gets 

= -^g^r''nxX'J^''Hq) (A13) 

where 
in which 

Q{q, a) = a(l - a)q^ - m^. (A15) 

Now, we arc in position to examine the imaginary parts of the functions 11^^^^,°'' and 11^^°''. We First write down 
the imaginary parts of the functions j'^' and J^'^\ 

ImJ^^^ = r^sm{xm^) (A16) 

(47r)2 7o ' 

and 

ImJ^''^ = — ^ da ^ sin{a;[a(l - a)q'' - m']}. (A17) 

(^i") Jo Jo 2x 

For the integral over a;, obviously, the major contribution arises from the integrand at the neighborhood of the origin. 
Therefore 

7mj(^) > 0. (A.18) 

As for the ImJ^"^^ , the integral over a may be estimated by taking the mean value | of the variable a in the integrand. 
Thus, 

It is well known that in the physical region, 

> Am^ (A20) 
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where = ^rr? is the starting point of a cut which is the solution of the following Landau equations'^^l 

Ai()fe2-m2) = 0, 

A2[((Z-fc)' -m2] =0, (A21) 
Aifc^ - A2(g - A:)^ = 0. 

In view of Eq. (A20), we may conclude 

/mj(2) > 0. (A22) 

The results in Eqs. (A18) and (A22) straightforwardly lead to that the imaginary parts of the vS^^^q) and '^^S^yf'iyq) 
have the same sign as we see from Eq. (A2) and (A13). 

At last , we would like to mention the imaginary part of the loop in Fig. (7c). The expression of the loop is 

ni'if '(g) = -9T''f'"'jtl («) (A23) 

where 

Completely following the procedure formulated in Eqs. (A8)-(A12), one may derive 

•^AA'(9) = J^l/^^f^^i^axx' - a{l - a)g.gv]e»'5(«.«). (A25) 
Neglecting the second term containing ^a^A' and then substituting the above expression into Eq. (A23), we can write 

nfi"' il) = -9^r''f"'9xx' J^') [q] ( A26) 

where 




(A27) 



Clearly, the sign of the imaginary part 

ImJ^''\q) = --—^ da ^ sm{x[a{l - a)q^ - m'']} (A28) 

is opposite to the ImJ^'^\q) shown in Eq. (A17). Therefore, the imaginary parts of the \l^^j^?^{q) and the \l^^y^{q) 
have opposite signs. 
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[1] FIGURE CAPTION 

Fig. (1): The tree diagram for fermion and antifermion scattering. 
Fig. (2): The tree diagram for two-gauge boson scattering. 

Fig. (3): The fourth-order diagrams for two-gauge boson scattering with only the gauge boson intermediate states. 
Fig. (4): The fourtii-ordcr diagrams for two-gauge boson scattering with only the ghost particle intermediate states. 
Fig. (5): The tree diagrams used to give all the diagrams in Fig. 3 through folding them with their conjugates 
Fig. (6): The tree diagrams used to give all the diagrams in Fig. 4 through folding them with their conjugates. 
Fig. (7a): The one-loop diagram formed by four-line gauge boson vertex. 
Fig. (7b): The one-loop diagram formed by three-lino gauge boson vertex. 
Fig. (7c): The one- loop diagram formed by ghost intermediate states. 
Fig. (8): The fourth-order diagrams for fermion-antifermion scattering. 

Fig. (9): The tree diagrams used to give all the diagrams in Fig. 8 through folding them with their conjugates. 
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